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Newtonian limit of Extended Theories of Gravity (in particular, higher- 
order and scalar-tensor theories) is theoretically discussed taking into account 
recent observational and experimental results. Extended Theories of Gravity 
have become a sort of paradigm in the study of gravitational interaction since 
several motivations push for enlarging the traditional scheme of Einstein Gen- 
eral Relativity. Such issues come, essentially, from cosmology and quantum 
field theory. In the first case, it is well known that higher-derivative theories 
and scalar-tensor theories furnish inflationary cosmological solutions capa- 
ble, in principle, of solving the shortcomings of The Standard Cosmological 
Model. Besides, they have relevant features also from the quantum cosmol- 
ogy viewpoint. In the second case, every unification scheme as Superstrings, 
Supergravity or Grand Unified Theories, takes into account effective actions 
where nonminimal couplings to the geometry or higher-order terms in the 
curvature invariants come out. Such contributions are due to one-loop or 
higher-loop corrections in the high-curvature regimes near the full (not yet 
available) quantum gravity regime. In the weak-limit approximation, all these 
classes of theories should be expected to reproduce the Einstein General Rel- 
ativity which, in any case, is experimentally tested only in this limit. This 
fact is matter of debate since several relativistic theories do not reproduce 
the Einstein results in the Newtonian approximation. For example, some au- 
thors claim for explaining the flat rotation curves of galaxies by Yukawa-like 
corrections to the Newtonian potential arising in this context. Others have 
shown that a conformal theory of gravity could explain the missing matter 
problem. Moreover, indications of an apparent, anomalous, long-range accel- 
eration revealed from the data analysis of Pioneer 10/11, Galileo, and Ulysses 
spacecrafts could be framed in a general theoretical scheme by taking into 
account Yukawa-like or higher-order corrections to the Newtonian potential. 
Such data will be very likely capable of confirming or ruling out the physical 
consistency of Extended Theories of Gravity. 
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I. INTRODUCTION 



Einstein's General Relativity (GR) can be considered one of the major scientific achieve- 
ment of last century. For the first time, a comprehensive theory of spacetime, gravity and 
matter has been formulated giving rise to a new conception of the Universe. However, in the 
last thirty years, several shortcomings came out in the Einstein scheme and people began to 
investigate if GR is the only and fundamental theory capable of explaining completely gravi- 
tational interaction. Such issues come, essentially, from cosmology and quantum field theory. 
In the first case, the presence of the Big Bang singularity, flatness and horizon problems 
[1] led to the statement that Standard Cosmological Model [2], based on GR and Standard 
Model of particle physics, is inadequate to describe the Universe in extreme regimes. On 
the other hand, GR is a classical theory which does not work as a fundamental theory, if 
one wants to achieve a full quantum description of spacetime (and then of gravity). Due 
to this facts and, first of all, to the lack of a definitive quantum gravity theory, alternative 
theories of gravity have been pursued in order to attempt, at least, a semi-classical scheme 
where GR and its positive results could be recovered. 

One of the most fruitful approach has been that of Extended Theories of Gravity (ETG) 
which have become a sort of paradigm in the study of gravitational interaction based on the 
enlargement and the correction of the traditional Einstein scheme. The paradigm consists 
in adding higher-order curvature invariants and minimally or non-minimally coupled scalar 
fields into dynamics which come out from quantum terms in the effective action of gravity. 
Other motivations come from the issue of a complete recovering of Mach principle [3] which 
leads to assume a varying gravitational coupling. All these approaches are not the '\full 
quantum gravit]/' but are needed as working schemes toward it. In any case, they are going 
to furnish consistent and physically reliable results. 

For example, it is well known that higher derivative theories [4] and scalar tensor theories 
[5] give rise to inflationary cosmological solutions capable, in principle, of solving the short- 
comings of Standard Cosmological Model. Besides, they have relevant features also from 
the quantum cosmology point of view since they give interesting solutions to the problem 
of initial, at least in the restricted context of minisuperspaces [6] . 

Furthermore, every unification scheme as Superstrings, Supergravity or Grand Unified 
Theories, takes into account effective actions where nonminimal couplings to the geometry 
or higher-order terms in the curvature invariants come out. Such contributions are due 
to one-loop or higher-loop corrections in the high-curvature regimes near the full (not yet 
available) quantum gravity regime [7]. However, in the weak-limit approximation, all these 
classes of theories should be expected to reproduce GR which, in any case, is experimentally 
tested only in this limit [8]. 

This fact is matter of debate since several relativistic theories (in particular ETG) do 
not reproduce exactly Einstein results in the Newtonian approximation but, in some sense, 
generalize them. In fact, as it was firstly noticed by Stelle [9], a i?^-theory gives rise to 
Yukawa-like corrections to the Newtonian potential which could have interesting physical 
consequences. For example, some authors claim for explaining the flat rotation curves of 
galaxies by using such terms [10]. Others [11] have shown that a conformal theory of gravity 
is nothing else but a fourth order theory containing such terms in the Newtonian limit and, 
by invoking these results, it could be possible to explain the missing matter problem without 
dark matter. 

Besides, indications of an apparent, anomalous, long-range acceleration revealed from the 
data analysis of Pioneer 10/11, Galileo, and Ulysses spacecrafts could be framed in a general 
theoretical scheme by taking into account Yukawa-like or higher order corrections to the 
Newtonian potential [12]. 

In general, any relativistic theory of gravitation can yield corrections to the Newton po- 
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tential (see for example [13]) which, in the post-Newtonian (PPN) formahsm, could furnish 
tests for the same theory [8]. 

Furthermore the newborn gravitational lensing astronomy [14] is giving rise to additional 
tests of gravity over small, large, and very large scales which very soon will provide direct 
measurements for the variation of Newton coupling Gn [15], the potential of galaxies, clus- 
ters of galaxies [16] and several other features of gravitating systems. Such data will be very 
likely capable of confirming or ruling out the physical consistency of GR or of any ETC 

In this review, we want to discuss the Newtonian limit of ETG, trying to establish a link 
to GR by putting in evidence the common features. In particular, we show that, in the 
Newtonian limit, such theories can have features which "generalize" the Newton potential 
and could explain several observational results. For example, a quadratic correction to the 
Newtonian potential strictly depends on the presence of a scalar-field potential which acts, 
in the low energy limit, as a cosmological constant. 

The layout of the review is the following. In Sec. II, we discuss the foundation of any metric 
theory of gravity: the Equivalence Principle in all its formulations. After, we develop the 
Parametrized Post Newtonian (PPN) formalism which constitutes the basis to compare every 
theory of gravity. Moreover, we outline the main features of GR discussing, in particular, the 
weak field approximation, the linearized Einstein equations and the gravitational radiation. 
Finally, a discussion of classical tests of GR is pursued considering the current experimental 
limits of the PPN parameters. 

Sec. Ill is devoted to the ETG which are classified according to their main features. We 
start with the minimal extension of GR which is the inclusion of a cosmological constant 
in Einstein equations. After, we discuss scalar-tensor and higher-order theories of gravity 
deriving, in every case, the field equations. 

In Sec. IV, the Newtonian limit of ETG is discussed considering, in particular, fourth-order 
and scalar-tensor theories. As an application, we develop the approach for String-Dilaton 
gravity deriving its weak field limit. 

The discussion of the results and the conclusions are drawn in Sec.V. 

II. THE FOUNDATION OF METRIC THEORIES: THE EQUIVALENCE 

PRINCIPLE 

Equivalence principle (EP) is the foundation of every metric theory of gravity [8]. 

The first formulation of EP comes out from the theory of gravitation studied by Galileo 
and Newton; it is called the Weak Equivalence Principle (WEP) and it states that the "in- 
ertial mass" m and the "gravitational mass" M of any object are equivalent. In Newtonian 
physics, the "inertial mass" m is a coefficient which appears in the second Newton law: 
F = ma where F is the force exerted on a mass to with acceleration a; in Special Relativity 
(SR) (without gravitation) the "inertial mass" of a body appears to be proportional to the 
rest energy of the body: E = mc^. Considering the Newtonian gravitational attraction, 
one introduces the "gravitational mass" M: the gravitational attraction force between two 
bodies of "gravitational mass" M and M' is F = GnMM' /r'^ where Gn is the Newtonian 
gravitational constant and r the distance between the two bodies. Various experiments [17] 
demonstrate that to = M. The present accuracy of this relation is of the order of 10~^^; 
spatial projects are currently designed to achieve precision of 10~^^ [18] and 10~^* [19]. 

The WEP statement implies that it is impossible to distinguish between the effects of 
a gravitational field from those experienced in an uniformly accelerated frame, using the 
simple observation of the free-falling particles behavior. The WEP can be formulated again 
in the following statement [8]: 
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If an uncharged test body is placed at an initial event in spacetime and given an 
initial velocity there, then its subsequent trajectory will be independent of its 
internal structure and composition. 

A generalization of WEP claims that SR is only locally valid. It has been achieved by 
Einstein after the formulation of SR theory where the concept of mass looses some of its 
uniqueness: the mass is reduced to a manifestation of energy and momentum. According 
to Einstein, it is impossible to distinguish between uniform acceleration and an external 
gravitational field, not only for free-falling particles but whatever is the experiment. This 
equivalence principle has been called Einstein Equivalence Principle (EEP); its main state- 
ments are the following [8]: 

a) WEP is valid; 

b) the outcome of any local non-gravitational test experiment is independent 

of velocity of free-falling apparatus; 

c) the outcome of any local non-gravitational test experiment is independent 
of where and when in the Universe it is performed. 

One defines as "local non-gravitational experiment" an experiment performed in a small- 
size * freely falling laboratory. From the EEP, one gets that the gravitational interaction 
depends on the curvature of spacetime, i.e. the postulates of any metric theory of gravity 
have to be satisfied [8] : 

a) spacetime is endowed with a metric g^^; 

b) the world lines of test bodies are geodesies of the metric; 

c) in local freely falling frames, called local Lorentz frames, the non-gravitational laws of 
physics are those of SR. 

One of the predictions of this principle is the gravitational red-shift, verified experimen- 
tally by Pound and Rebka in 1960 [20]. 

It is worth stressing that gravitational interactions are specifically excluded from WEP 
and EEP. In order to classify alternative theories of gravity, the Gravitational Weak Equiva- 
lence Principle (GWEP) and the Strong Equivalence Principle (SEP) has to be introduced. 
The SEP extends the EEP by including all the laws of physics in its terms [8] : 

a) WEP is valid for self-gravitating bodies as well as for test bodies (GWEP); 

b) the outcome of any local test experiment is independent of the velocity of 

the free-falling apparatus; 

c) the outcome of any local test experiment is independent of where and when 
in the Universe it is performed. 

Therefore, the SEP contains the EEP, when gravitational forces are ignored. Many au- 
thors claim that the only theory coherent with the SEP is General Relativity (GR). 



*In order to avoid the inhomogeneities. 
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A. The Parametrized Post Newtonian (PPN) limit 



GR is not the only theory of gravitation and, up to now, at least twonty-five alternative 
theories of gravity have been investigated from the 60's, considering the spacetime to be 
"special relativistic" at a background level and treating gravitation as a Lorentz-invariant 
field on the background. 

Two different classes of experiments have been studied: the first ones testing the foun- 
dations of gravitation theory ~ among them the EP - the second one testing the metric 
theories of gravity where spacetime is endowed with a metric tensor and where the EEP is 
valid. However, for several fundamental reasons which we are going to discuss, extra fields 
might be necessary to describe the gravitation, e.g. scalar fields or higher-order corrections 
in curvature invariants. 

Two sets of equations can be distinguished [17]. The first ones couple the gravitational 
fields to the non gravitational contents of the Universe, i.e. the matter distribution, the elec- 
tromagnetic fields, etc... The second set of equations gives the evolution of non-gravitational 
fields. Within the framework of metric theories, these laws depend only on the metric: this is 
a consequence of the EEP and the so-called "minimal coupling" . In most theories, including 
GR, the second set of equations is derived from the first one, which is the only fundamental 
one; however, in many situations, the two sets are decoupled. 

The gravitational field studied in these approaches (without cosmological considerations) 
is mainly due to the Sun and the Eddington-Robertson expansion gives the corresponding 
metric. Assuming spherical symmetry and a static gravitational field, one can prove that 
there exists a coordinate system (called isotropic) such as 



dt being the proper time between two neighboring events. 

The Newtonian gravitational field does not exceed GnMq/Rqc^ ^ 2 x 10~^, where c 
is the speed of light, Mq is the mass of the Sun and Rq its radius. The metric is quasi- 
Minkowskian, A{r) and B{r) are dimensionless functions which depend only on Gn, M, 
c and r. Indeed, the only pure number that can be built with these four quantities is 
GnM/tc^. The Eddington-Robertson metric is a Taylor expansion of A and B which gives 



The coefficients a, /3, 7 are called the post Newtonian parameters and their values depend 
on the considered theory of gravity: in GR, one has a = /? = 7 = 1. The post Newtonian 
parameters can be measured not only in the Solar System but also in relativistic binary 
neutron stars such as PSR 1913 -|- 16. 

A generalization of the previous formalism is known as Parametrized Post-Newtonian 
(PPN) formalism. Comparing metric theories of gravity among them and with experi- 
mental results becomes particularly easy if the PPN approximation is used. The following 
requirements are needed: 

• particles are moving slowly with respect to the speed of light; 

• gravitational field is weak and considered as a perturbation of the flat spacetime; 

• gravitational field is also static, i.e. it does not change with time. 



ds^ = B{r) d^ - A{r) - (rfr^ sin^ e dcj)'^ + dO"^) . 



(1) 
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The PPN limit of metric theories of gravity is characterized by a set of 10 real-valued 

parameters; each metric theory of gravity corresponds to particular values of PPN parame- 
ters. The PPN framework has been used first for the analysis of Solar System gravitational 
experiments, then for the definition and the analysis of new tests of gravitation theory and 
finally for the analysis and the classification of alternative metric theories of gravity. 

By the middle 1970's, the Solar System was no more considered as the unique testing 
ground of gravitation theories. Many alternative theories of gravity agree with GR in the 
Post-Newtonian limit and thus with Solar System experiments: nevertheless, they do not 
agree for other predictions (such as cosmology, neutron stars, black holes and gravitational 
radiation) for which the post-Newtonian limit is not adequate. In addition, the possibility 
that experimental probes, such as gravitational radiation detectors, woiild be available in 
the future to perform extra-Solar System tests led to the abandon of the Solar System as 
the only arena to test gravitation theories. 

The study of the binary pulsar PSR 1913 + 16, discovered by R. Hulse and J. Taylor 
[21], showed that this system combines large post- Newtonian gravitational effects, highly 
relativistic gravitational fields (associated with the pulsar) with the evidence of an emission 
of gravitational radiation (by the binary system itself). Relativistic gravitational effects 
allowed one to do accurate measurements of astrophysical parameters of the system, such 
as the mass of the two neutron stars. The measurement of the orbital period rate of change 
agreed with the prediction of the gravitational waves (GW) emission in the GR framework, 
in contradiction with the predictions of most alternative theories, even those with PPN 
limits identical to GR. However, the evidence was not conclusive to rule out other theories 
since several shortcomings remain, up to now, unexplained. 



B. Basics of General Relativity 

GR is currently the most successful relativistic theory of gravitation, introduced by A. 

Einstein [22] in 1915 to include in a natural way the EP and SR theory (1905). Einstein's 
derivation of GR has been driven by theoretical criteria of elegance and simplicity, not in- 
volving neither experimental nor observational results. Nevertheless, he had soon to compare 
his theory with experimental tests. Three important ones confirmed GR: i) the anomalous 
perihelion of Mercury (1915) [8]; ii) the defiection of light by the Sun (1919) [23]; Hi) the 
gravitational red-shift of light - measured by Pound- Rebka experiment (1960) [20]. GR is 
a metric theory where the metric tensor is the only gravitational field [2,24-26]. 
The Einstein tensor G^i^ is defined as 

G^v = Ri^v — -^dn'^^ i^) 

with R^i, being the Ricci tensor computed by contracting the Riemann tensor i?^,^ — R'^^^^ 
and R = R^^ = g^^Rfiv the curvature scalar. The Einstein tensor is symmetric thanks to the 
features of Ricci and metric tensors. It verifies the Bianchi identity V^G^v = 0. The second 
term in the r.h.s. of Eq. (3) is divergence-free as the stress-energy tensor T^j^ because of the 
conservation of energy-momentum. 

Einstein equations present, in a formal way, the fundamental relation characterizing the 
spacetime description of GR: matter curves spacetime around it. 

Such a relation connects two 4-dimensional symmetric tensors of rank 2, through a propor- 
tionality relation: the Einstein tensor G^^, defined above as a contraction of the Riemann 
tensor Ra/s-yS, describing the curvature of a generic manifold, and the stress-energy ten- 
sor T^i, which contains densities and energy-momentum fluxes. This tensor depends on all 
possible types of energy (kinetic, potential, pressure, etc...) and also on gravitational field 
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energy. Therefore, Einstein equations are non-linear since gravity terms appear on both 
sides (source and effect). 

Einstein equations can be obtained in two different ways. The first one is the natural 
covariant generalization of Poisson equation for the Newtonian gravitational potential, i.e. 



G^. = X T^. with X = — ^ • (4) 



The proportionality constant x is computed by assuming the convergence to the Newto- 
nian theory of gravity in the limit of weak, time-independent gravitational fields and slow 
motions. 

Einstein equations arc non-linear second-order differential equations; there are ten inde- 
pendent equations since both sides are two-index symmetric tensors (the solution consists in 
computing the ten unknown functions of the metric components); nevertheless, the Bianchi 
identities give four additional constraints on R^j^i, . so there arc only six independent relations 
in Eq. (4). It is not possible to solve these equations in general, but even particular solutions, 
like the ones in vacuum, are difficult to compute without some simplifying assumptions. 

The second way to derive Einstein equations comes from the variation of the action prin- 
ciple where the curvature scalar invariant R is considered. The GR action is the Hilbert- 
Einstein action [25] 

^ Jd^Xy/^[R+Cm] (5) 



16wG 

linear in the Ricci curvature scalar R, minimally coupled, being the gravitational coupling 
the constant On, where g is the metric determinant and £m is the perfect-fluid matter 
Lagrangian density. 

Assuming spherical symmetry - Eq. (1) - there are standard coordinates such that 

ds'^ = e^dt'^-e^dr'^-r'^ism^edcp'^ + de^) (6) 

where u and A are two functions of t and r only. 
The solution of the Einstein equations in vacuum is 

e'' = e~^f{t) and = (7) 

where f{f,) is an arbitrary positive function of f, and a is a constant. The metric is time - 
independent: this is the Birklioff's theorem. For instance, it states that the explosion 
of a supernova does not produce any time dependent modification of spacetime if it is a 
spherically symmetric explosion, i.e. there is no GW production ^. 

Let us assume that the source of gravitational field is a mass located at the origin r = 0. 
In order to recover the Newtonian limit far away from the origin, one must assume a/r = 
2(j) = — 2(Gjv-M)/(rc^). The corresponding solution is called the "Schwarzschild solution" 
(in standard coordinates) 

ds^ = B{r) dt"--^- (sin^ 6 dcj)'' + d^^) ^jth B{r) = 1 - i^iL^ (8) 
B[r) r 

where M is the mass of the central body. 

Generally, the Schwarzschild radius of the central mass, Rs = 2Gn M/c?, is much smaller 



This is a significant limitation for the GW production as large asymmetries are rather rare. 



6 



than the physical radius i?o of the body itself. Therefore the solution (8) is only valid in 

a region where r > Rq > Rs- However, one can also consider a compact object with 
Ro < Rs', in that case the surface r = Rs separates the spacetime into two regions. Such a 
compact object is called a Schwarzschild black hole. In the exterior region, where r > Rs, 
the coordinates t and r arc respectively time like and space like (gti = B > and g^r = 
— 1/B < 0). When a light source approaches the surface r = Rs, the Einstein gravitational 
shift becomes so huge that the frequency of the light becomes zero everywhere; the photons 
loose their energy and cannot be seen anymore. The surface r = Rs is the last surface that 
can be seen from outside and it is thus called the "horizon" . In the interior region, where 
r < Rs one has gu = B <0 and grr = > 0; it means that t is a space-like coordinate 

while r is a time-like coordinate. 

1. The weak field approximation and the linearized Einstein equations 

Replacing the Newtonian limit by a less restrictive hypothesis leads to the weak field 

approximation: practically, the field is still weak, but it is allowed to change in time and 
there is no more restriction on the test particles motion. New phenomena are associated 
with this hypothesis like the emission of gravitational radiation and the deflection of light. 
This framework allows one to split the metric into two parts: the flat Minkowski metric 
ViJ.i' = diag(c, —1, —1, —1) plus a perturbative term h/j,,^, assumed to be small. This linearized 
version of GR describes the propagation of a symmetric tensor h^i, on a flat background 
spacetime. So, the metric reads 

g^iy = tlnv + h^,^ with \h^,v\^l. (9) 

As h is small, one can neglect terms higher than the flrst order in the perturbation h^v; 
in particular, one can raise/lower indexes with ry^y and rj^" as the corrections are of higher 
order in the perturbation 

gt^u ^ _ f^^u ^.^jj f^^.^ ^ V^^V'""hpa . (10) 

The aim is to find the equations of motion to which the perturbations /i^,^ obey by 
investigating the Einstein equations to the first order. Inserting the new metric (9) in the 
Einstein tensor, we obtain 

Gfiu = ^ {dadt,h'^ + dadij,h'l - d^duh- hij,t, - rjnt.dpdah'"^ + 'iji^uh) (11) 

where h = rj^^^h^y = h'^ is the trace of the perturbation and = ^dtt — dxx — dyy — dzz is the 
d'Alembertian of the flat spacetime, using from now on (unless otherwise stated) geometri- 
cal units for which c = 1. 

The stress-energy tensor is computed at the 0-order in hixiy' the energy and the mo- 
mentum have to be small too, according to the weak field approximation and the lowest 
non-vanishing order in is of the same order of magnitude as the perturbation. Therefore, 
the conservation law becomes d'^T^i, = 0. 

2. Gravitational Waves 

GW are weak ripples in the curvature of spacetime, produced by the motions of matter. 
They propagate at the speed of light. The linearized Einstein equations allow wave solutions, 
in a way similar to electromagnetism. These GW are transverse to the propagation direction 
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and show two independent polarization states. 

The new metric (9) docs not fix the spacetimc frame completely; two possible gauges 
can be applied in addition to simplify the Einstein equations. Using the Lorentz gauge 
d^h'^ — \d\h = 0, the Einstein equations (4) are linearized and can be written as 

Oh^,,-^r]^,ah = . (12) 

The trace-reversed perturbation is defined as 

h^,y = hf,i,-^ riu,y h . (13) 

One can choose a frame in which the harmonic gauge condition, d^h'^ = 0, is verified. 
Then, the Einstein field equations (4) become 

□V = ^^V- (14) 

and, in vacuum, one has simply 

□V = 0- (15) 

These equations are similar to the electromagnetism field equations and one can use the 
same method to solve them. Indeed, looking at the Einstein equations in vacuum Eq. (15) 
-, one can note that they are in the form of a wave equation for /i^y - the D'Alembertian 
reduces to the form □ = ^du — V^. Therefore, in the absence of matter, one looks for plane 
waves solutions 

V = C^^e''^'-' (16) 

where C^,y is a constant and symmetric tensor of rank 2 and A; is a time-like vector, the 
wave vector. The plane wave in Eq. (16) is a solution of the linearized equations in vacuum 
- Eq. (15) - if the wave vector k is null, i.e. satisfies the condition Wk^ = and shows that 
GW propagate to the speed of light. 

The four conditions of the harmonic gauge fc^C^'' = lead to six independent components 
for the symmetric tensor C^v As there are still some unused degrees of freedom, one can 
make another gauge choice on the tensor C^^: 

C/^ = Traceless; (17) 
C°'^ = Transverse. (18) 

One has, in this way, the so called Transverse- Traceless (TT) gauge. These five relations 
give four new constraints on in addition to the harmonic gauge condition; therefore, 
only two independent components remain in C^" . As the wave is traceless, one can check 
from Eq. (13) that /i^J = /iJJ. 

Therefore, the general form of the symmetric tensor C^v is finally 
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(19) 



Let us define Cu — h+ and C12 — hx] GW appears to have two polarized states - 
and "x" - which modify the spacetime curvature in different ways. In tensorial form, one 
can write 
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h = [h+ (el 8) el - e"^ e"^) + 2 /ix (el ea)] e''^(*-^/'=) 



(20) 



Being ^ = ($^,^^,^^) the separation between two free particles and taking into account 
the geodesic deviation [24] which describes the evolution of two free- falling particles, if the 

GW propagates in the direction x^, only and arc involved in its passage. Assuming a 
polarized GW, the integration of the geodesic deviation equation gives: 

• Polarization "+" {h^ = 0) : 

f eit)\ _ f i + ^h+e"'''=' 

• Polarization "x" {h+ = 0) : 




Let us consider now a test-mass ring (a massless and free-falling set of particles) inter- 
acting with GW, lying in an plane orthogonal to the direction of the wave propagation. Its 
oscillations depend on the GW polarization. 

After having found a solution to Einstein equations in vacuum, let us solve now Eq. (14) 
with a non-zero source term. The solution is computed using the retarded Green function 

V(i,^) = ^ / —^T^,itr,y)d^y (23) 

^ Jye Source I V \ 

with \ x — y\ = \/Sij{x^ — y*){x^ — y^) (Euclidean distance) and tr = t—\x — y \ /c (retarded 
time). 

Let us consider an isolated source with a density p and a characteristic dimension 8R, 
located at a distance R from the observation point x. One assumes 5R <^ R so, in particular, 
one gets \ x — y\ !^ R and one can move this constant term outside the integral in Eq. (23). 
As the stress-energy tensor verifies the conservation of energy df^T'^" = 0, the harmonic 
gauge condition d^h'^ = is also verified. Moreover, the radiation is mostly emitted at 

5R 1 

frequencies cj/27r, so that — <C — . Then, it is possible to demonstrate that only the 

C (jj 

spatial coordinates of tensor /i^,^ are difi^erent from zero. 

The quadrupole momentum tensor of the source energy density is defined as 

quit) = [ yiyjToo{t,y}(fy with Too « pc^ . (24) 

ye source 

The metric perturbation is given by 

hij {t, x) = Qij (tr) . (25) 

So, GW, generated by an isolated non-relativistic object, are proportional to the second 
derivative of the quadrupolar momenta of the energy density. Eq. (25) shows that the 
metric perturbation amplitude h varies as the inverse of distance to the source R; a faster 
decreasing with the distance, e.g. 1/R?, would make vain the hope of any GW detection. 
Fortunately, GW detectors arc sensitive to the amplitude of the signal. 

The energy emitted by gravitational radiation is difiicult to define. A way to overcome 
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this difficulty is to define the stress-energy tensor by developing the metric g^v and the 
Einstein tensor G^i, at the second order: 

g^iu = -n^v + Ku + h^^} , (26) 

= G« [v + h^'^] + G^i) [v + h]. (27) 

Einstein equations in the vacuum = can be written in the form 

Gi'}[v + h^'^] = (28) 



with the definition 



c4 



t^'^ = + (29) 

» TT Cr ^ 



N 



The Bianchi identity says that d^J.^^ = 0, therefore t^^ can be considered as the stress- 
energy tensor of a gravitational field - yet, it is only a pseudo-tensor [25]. One can compute 
the energy density too by averaging over many cycles (because the energy is not localized) 
the GW energy: 



C , ,• 2 -2 



too = +/ix ) . (30) 



Then, the source luminosity is L = r / ctoodfl. Introducing the reduced 

J Sphere of radius r 



quadrupolar momenta 



Qij = Qij - ^Sij6'''-qki (31) 



one obtains the Einstein quadrupole formula 

It is worth noticing that a GW emission requires a variation of the quadrupolar mo- 
mentum, as shown already in Eq. (25). This typical feature of gravitational radiation, 
added to the weakness of the coupling constant between gravitation and matter - here 
G/c^ w 10~^^ - explains why GW amplitudes are so small compared to those pro- 

duced by electromagnetic radiation. In conclusion, the gravitational radiation is quadrupolar 
and a symmetric spherical body does not emit GW because its reduced quadrupolar mo- 
menta are zero. 

The corresponding quantum field is the graviton with zero mass and spin 2. Within the 
framework of more general theories, the gravitation can be described as a combination of 
states with a spin 2 and spin zero or as a particle with spin maximum 2. One can also 
imagine that the mass of graviton is not zero or that the state mass spectrum is complex. 
Presently, there is no observational reason to doubt that the present observational bounds 
on the mass of the graviton are severe. 

Discovering GW would open "a new window onto the Universe" [27]. It is clear that being 
sensitive to an additional radiation would lead to major discoveries like when the Universe 
became observed through radio, X or gamma waves. Then, it would allow physicists to test 
GR in the strong field regime, to check the gravity velocity (assumed to be the speed of 
light in the Einstein theory) or to verify that GW only change distances perpendicular to 
their direction of propagation. 
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Alternative relativistic theories of gravity also predict the existence of GW. However, 

many essential features of the radiation arc different: the number of polarization states, the 
propagation speed, or the efSciency of wave generation. We shall not take into account GW 
any further in this review, but it is worth stressing that their detection is of fundamental 
interest since their features and parameters, through the PPN limit or in strong field regime, 
allow a discrimination among relativistic theories of gravity [27] . 



Historically, the gravitational red - shift experiment, the deflection of light and the anoma- 
lous perihelion shift of Mercury are considered as the three classical tests of GR [8]: indeed, 
they were among the first observable effects of GR to be computed by Einstein. The grav- 
itational red-shift experiment can be considered as a test of the EEP. upon which GR and 
other metric theories of gravity are founded. Moreover, a further test has been introduced: 
the time delay of light due to the presence of gravitational field. 

The deflection of light by the Sun is described, in physical units, by the deflection angle 



where is the angle between the unperturbed paths of the photons from the source, the 
Sun, and the reference source and D is the distance between the center of mass and the 
point of closest approach of the unperturbed ray. For 7 = 1 one finds the results of GR. 

The explanation of the anomalous perihelion shift of Mercury orbit was the subject of a 
controversy to establish if this shift was a confirmation or a refutation of GR because of the 
apparent existence of a Solar quadrupole momentum that could contribute a fraction of the 
observed perihelion shift. 

The time delay of light depends on the strength created by the gravitational field of a 
massive body located near the photon path: because of this effect, a light signal will take a 
longer time to cross a given distance than it would if the Newtonian theory was valid. The 
time delay is 



where M is the mass of the body, Rq the Earth radius, Rq the Solar radius and Rp the 
distance of a planet or a spacecraft. 

Beside experimental tests, it is worth mentioning some theories about tests: kinematic 
frameworks [28-31] to prove the Local Lorentz Invariance (LLI), dynamical frameworks 
[32,33] to test WEP, LLI and Local Position Invariance (LPI), unification theory [34] moti- 
vated to test WEP, LLI and LPI. 



C. Discussion of Classical Tests of General Relativity 




(33) 




(34) 
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Tab. I shows the PPN parameters and their significance in GR, in semi-conservative the- 
ories and in scalar-tensor theories [8]. 
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A direct test of the WEP is the Ebtvos experiment, i.e. the eomparison of the rest accel- 
eration of two laboratory-sized bodies of different composition in an external gravitational 
field. Indeed, the accelerations of different bodies are not equal if the WEP is not valid. A 
measurement or limit on the relative difference in acceleration yields to a quantity called 
^^Eotvos ratio" given by 

2^1^ (35) 
\ai+a2\ ^ ' 

where a\ and a-i are the accelerations of the two bodies. 

Experimental limits on r\ place limits on the WEP-violation parameter. Several high- 
precision Ebtvbs-ty\tc experiments have been performed and limits on from different ex- 
periments are shown in Tab. 11 [8]. 



Experiment 


Method 


Limit on |?7| 


Newton 
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10"^ 
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Boulder 


flotation on water 


4 X 10-" 


Orbital 


free fall in orbit 


IQ-i^ 10-1* not yet performed 
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Very few experiments have been used to make quantitative tests of the LLI in the same 

way that Eotvos experiments have been used to test the WEP. There is one experiment that 
can be interpreted as a "clean" test of the LLI: the Hughes-Drever experiment, performed 
independently in 1959 -=-60. If the LLI is violated, then there must be a preferred rest frame, 
presumably the one associated to the Cosmic Microwave Background (CMB), in which the 
local laws of physics take their special form. Deviations from this form would depend on the 
velocity of the laboratory relative to the preferred frame. New experimental tests will be 
performed in the future, using atomic clocks, macroscopic resonators or laser intcrferomctry, 
for instance the LATOR mission [35] which aims to carry out a test of the curvature of the 
Solar System's gravity field with an accuracy of 10~^. The two principal tests of LPI 
are gravitational red-shift experiments testing the existence of spatial dependence on the 
outcomes of local experiments, and the measurements of the possible time variations of the 
fundamental non-gravitational constants. Some Solar System experiments testing the SEP 
are based on the Lunar-Laser- Ranging (LLR) Eotvos experiment which takes into account a 
potential Nordtvedt effect (not yet observed), i.e. the violation of the GWEP. The breaking 
of the WEP for massive self-gravitating bodies (GWEP), which several metric theories 
predict, has a variety of observable consequences: the most important one is a polarization 
of the Moon orbit about the Earth. Because the Moon self-gravitational energy is smaller 
than the Earth one, the Nordtvedt effect causes that Earth and Moon fall toward the Sun, 
with slightly different acceleration. The "Nordtvedt" parameter rj and its limits are shown 
in Tab.III Its expression isr/ = 4/3-7-3 - (10/3)^ - ai (2/3)a2 - (2/3)Ci - (1/3)C2- In 
GR, rj = 0. 

Tests to establish preferred-frame and preferred-location effects can be divided into two 
main categories: geophysical tests and orbital tests. 

• Geophysical tests: some metric theories of gravity predict preferred-location ef- 
fects considering the locally estimated gravitational constant Gl, measured using 
Cavendish like experiments. These effects are violations of SEP. Limits are shown 
in Tab.III. 

• Orbital tests: there are a number of observable effects of preferred- frame and preferred- 
location types in the orbital motions of bodies governed by the n-body equations 
of motion. The most important ones are the perihelion shifts of planets coming in 
addition to the "classical" shift. 
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The current PPN parameter limits are shown in Tab. III. 
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Fig.l shows some experimental tests performed during a period of eighty years about the 
/I + 7\ 

parameter — - — . It is possible to derive a bound on w > 20, 000 in scalar-tensor gravity 
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FIG. 1. Measurements of the coefRcient (1 + "/)/2 from light deflection and time delay measure- 
ments - the GR value is unity. Shapiro time-delay measurements using Viking spacecraft and Very 
Large Baseline Interferometry (VLBI) light deflection measurements are in agreement with GR to 
0.1%. Results from Cassini mission are in agreement at the level of 5 x 10~^. 



From a quick view of limits and data in the tables, it is clear that no definite answer 
can be given on which is the " true" theory of gravity. GR works very well in weak field 
limit but shortcomings arise as soon as quantum fields, cosmology or strong gravitational 
fields are considered. Moreover, no definite tests have been performed outside Solar System 
and the "facf that GR works at larger scales {e.g. galactic sizes) is nothing else but a 
pure assumption (issues as "missing" or "dark" matter could be read as "indications" that 
GR does not work and further ingredients are needed to explain dynamics). Due to these 
reasons, alternative theories of gravity have to be studied and, possibly, tested. 
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III. EXTENDED THEORIES OF GRAVITY 



Extended Theories of Gravity (ETG) are, by definition, all the approaches to the physics 
of gravitational interaction which enlarge the Einstein scheme through an addition of cor- 
rective terms into the Hilbert-Einstein action (5). There are several reasons to do so and 
the resulting field equations are modified leading to different physical predictions. 

From a conceptual point of view, there is a priori no reason to restrict the gravitational 
Lagrangian to a linear function of the Ricci scalar i?,, minimally coupled with matter [36]. 
Furthermore, the idea that there are no "exact" laws of physics but that the Lagrangians 
of physical interactions are "stochastic" functions - with the property that local gauge in- 
variances {i.e. conservation laws) are well approximated in the low energy limit and that 
physical constants can vary - has been taken into serious consideration - see Ref. [37] for 
instance. This scheme was adopted in order to deal with the quantization on curved space- 
times and the result was that the interactions among quantum scalar fields and background 
geometry or the gravitational self-interactions yield corrective terms in the Einstein-Hilbert 
Lagrangian [38]. Moreover, it has been realized that such corrective terms are inescapable 
if wc want to obtain the effective action of quantum gravity on scales closed to the Planck 
length [39]. Higher-order terms in curvature invariants (such as R^, R^^^R^i,, R^^"^ R^iuaff^ 
RUR, on RU^'R) or non-minimally coupled terms between scalar fields and geometry (such 
as ip'^R) have to be added in the effective Lagrangian of gravitational field when quantum 
corrections are considered. For instance, one can note that such terms occur in the effective 
Lagrangian of strings or in Kaluza-Klein theories, when the mechanism of dimensional re- 
duction is used [40]. 

Using a completely different point of view, these extended theories become interesting 
when one tries to incorporate the Mach principle in gravity: this principle states that the 
local incrtial frame is determined by some average of the motion of distant astronomical ob- 
jects [41], so that gravitational coupling can be scale-dependent and related to some scalar 
field. As a consequence, the concept of "inertia" and equivalence principle have to be re- 
vised. For example, the Brans-Dicke theory is a serious attempt to define an alternative 
theory to the Einstein gravity: it takes into account a variable Newton gravitational con- 
stant, whose dynamics is governed by a scalar field non-minimally coupled with geometry. 
In such a way, the Mach principle is better implemented [3,42,43]. 

Beside fundamental physics motivations, all these theories have acquired a huge interest 
in cosmology due to the fact that they "naturally" exhibit inflationary behaviours able to 
overcome the shortcomings of Standard Cosmological Model (based on GR). The related 
cosmological models seem very realistic and, several times, capable of matching with the 
observations [4,5]. 

Furthermore, it is possible to show that, via conformal transformations, the higher-order 
and non-minimally coupled terms always correspond to the Einstein gravity plus one or 
more than one minimally coupled scalar fields [44-48]. More precisely, higher-order terms 
always appear as a contribution of order two in the equations of motion. For example, a term 
like R? gives fourth order equations [49], R Di? gives sixth order equations [48,50], RO^R 
gives eighth order equations [51] and so on. By a conformal transformation, any 2nd-order 
term corresponds to a scalar field: for example, fourth-order gravity gives Einstein plus one 
scalar field, sixth order gravity gives Einstein plus two scalar fields and so on [48,52]. 
This feature results very interesting if we want to obtain multiple inflationary events since 
an early stage could select "very" large-scale structures (clusters of galaxies today) , while a 
late stage could select "small" large-scale structures (galaxies today) [50]. The philosophy 
is that each inflationary era is connected with the dynamics of a scalar fleld. Furthermore, 
these extended schemes naturally could solve the problem of " graceful exit" bypassing the 
shortcomings of former inflationary models [5,53]. 
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In this section, we want to discuss some of these theories, outUning their main features. 
In particular, we shall consider: 

• GR with a cosmological constant, the minimal extension discussed by Einstein himself; 

• higher-order-scalar-tensor theories, which are the general scheme embracing all the 
ETG; 

• scalar-tensor theories of gravity, discussing, in particular, the most famous example 
which is the Brans-Dicke theory; 

• fourth-order theory of gravity; 

• higher than fourth-order theories of gravity. 

As a general feature (apart the first case), we shall see that the Einstein field equations 
result modified in two senses: i) geometry can be nonminimally coupled to some scalar 
field, and/or ii) higher than second order derivative terms in the metric come out. 

A. General Relativity with a cosmological constant 

The determination of the cosmological constant [54] is one of the most challenging issue 
of modern physics since by fixing its value, one could contribute to select self-consistent 
models of fundamental physics and cosmology. Briefly, its measurement should provide the 
gravity vacuum state [55] , should allow to understand the mechanism which led to the today 
observed large scale structure from the early universe [1,56], and to predict what will be the 
fate of the whole universe ("no-hair conjecture" - see below) [57]. 

The presence of a cosmological constant A modifies the Einstein equations (4); this term 
has been first included by Einstein himself and tuned to obtain a static Universe model. 
Later he rejected this "fine tuning" as his "biggest blunder" after the discovery by Hubble 
of the expansion of the Universe by looking at red-shifted galaxy spectra. Ironically, recent 
measurements [58,59] seem to favor a non-zero value of A which could be the major contri- 
bution to the cosmological density. 

Assuming the existence of a cosmological constant A, Einstein equations (4) can be written 

as 

Gu.iy = Tfj,„ + Ag^^ (36) 

and the action in presence of a cosmological constant is simply 

5 = -^/d^x^[E+2A]. (37) 

On the right-hand side of the Einstein equations, this A-term is a stress-energy tensor 
contribution automatically conserved by definition of the metric connection. It can be 
considered as a vacuum energy density, a source of energy and momentum present also in 
absence of matter fields; this interpretation agrees with quantum field theory which predicts 
that the vacuum should have some sort of energy and momentum. However, this prediction 
is 50 orders of magnitude above the current measurements of A! 

Thanks to the recent measurements of distant supernovae la luminosity [58], a non-zero 
cosmological constant A is coming back in the footlights. Nevertheless, other possibilities 
are still investigated to explain the origin of this "dark energy" , like the quintessence. Let 
Qa be the energy density associated to the cosmological constant, normalized to the critical 
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density of the Universe. Taking into account the recent results of WMAP experiment [59], 
its value is 



simeq 0.73 ± 0.04 (38) 



with the hypothesis Ototai = 1, well confirmed by experimental data. 
Prom the Eq. (38) it is possible to deduce a value for A 

A = (1.29 ±0.15) 10-^2 f——^:^^^——) m-2 (39) 

where Hq is the Hubble constant predicted by WMAP experiment [59]. The ratio between 
the gravitational force per mass unit due to a mass M, located at a distance £), and therefore 
due to the cosmological constant is given by [60] 

GjvM ^ ,„2i Earth-Sun distance \V \ 

as the total gravitational force per unit mass is F = — — "~"q — • '^■'^^ cosmological 

constant A is not negligible for cosmological distances. 

The no-hair conjecture [57] claims that if there is a positive cosmological constant, models 
of expanding universe will approach a de Sitter behavior, i.e. if a cosmological constant is 
present, the universe will necessarily become homogeneous and isotropic whatever its initial 
conditions were. However, there is no general proof of such a conjecture and there are even 
theoretical counter-examples of initially expanding and then re-collapsing universes which 
never become de Sitter [61]. 

On the other hand, a simplified version of the conjecture can be proved [62] : 

All Bianchi cosmologies (except Bianchi IX), in presence of a positive cosmolog- 
ical constant, asymptotically approach the de Sitter behaviour. 

It is worthwhile noticing that the cosmological constant is a true constant (added by 
hands) and that the contracted Bianchi identity is not used in the proof, which is thus 
independent of the matter evolution. In order to extend the "no-hair conjecture" to ex- 
tended theories of gravitation, one has to introduce different sets of conditions (with respect 
to those given in [62]) since the cosmological constant is not given, a priori, but can be 
recovered from the dynamics of scalar fields (considering higher-order geometric terms in 
the gravitational Lagrangian also as a kind of scalar fields [46,52]). 

Such conditions must not use the energy conditions [63] , but they have to allow the intro- 
duction of a kind of effective cosmological constant A^ff, which asymptotically becomes the 
de Sitter constant [54,64]. This feature is due to the fact that, in an expanding universe, all 
the contributions to the energy density and to the Riccd tensor must decay as some power 
of the scale factor a{t), while the cosmological constant is the only term which does not 
decrease with time. Hence, in an expanding universe, A is the asymptotically dominant 
term in the Einstein equations giving rise to a de Sitter spacetime. Actually, the effective 
cosmological constant is time-dependent but, at the end, it has to coincide with the de Sitter 
one (i.e. the true gravitational vacuum state). Then, given any extended theory of gravity, 
it could be possible, in general, to define an effective time-varying cosmological constant 
which becomes the "true" cosmological constant if and only if the model asymptotically 
approaches de Sitter {i.e. if the no-hair conjecture is recovered asymptotically). This fact 
introduces constraints on the choice of gravitational coupling, of scalar field potential and of 
higher-order geometrical terms whose combinations are used as components of the efi^ective 
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stress-energy tensor. Terms like R'^, OR and (p can be all treated as "scalar fields" in the 
construction of Ae//, i.e. all of them give rise to extra-terms in the field equations which 
contribute to the construction of an effective stress-energy tensor T'lfJ^\ 

In this context, a fundamental issue is to select the classes of gravitational theories and 
the conditions which "naturally" allow to recover an effective time-dependent cosmological 
constant without considering special initial data. 

In the framework of extended gravity theories, the meaning of a time dependent "cosmo- 
logical constant" can be discussed [54]. As we said, the main request is that such effective 
cosmological constants evolve (at least asymptotically) toward the actual cosmological con- 
stant which means that the de Sitter behaviour has to be recovered. 



B. Higher-order-scalar-tensor theories: the general scheme 



If the gravitational Lagrangian is nonlinear in the Ricci scalar (and, in general, in the 
curvature invariants) its dynamics, {i.e. the Einstein equations), has an order higher than 
second; for this reason, such theories arc often called higher-order gravit,ational theories. 
Physically, as we said, they are interesting since higher-order terms in curvature invariants 
appear when one performs a one-loop renormalization of matter and gravitational fields in 
a curved background (sec for example Rcf. [38,65]). 

We can take into account the the most general class of higher-order-scalar-tensor theories 
in four dimensions. They can be generated by the action 



d^x- 



F{R,UR,U^R,..U^R,ct>)--, 



(41) 



where F is an unspecified function of curvature invariants and of a scalar field (p. The 
term C^n is, as above, the minimally coupled ordinary matter contribution. Physical units, 
SttG = c = h = 1, will be used, unless otherwise stated; e is a constant which specifies the 
kinetic term of the theory. 

The field equations are obtained by varying Eq. (41) with respect to the metric £f^j/. One 
gets 



T^-" + \ 9^" [F - GR) + {g^^g-^ - g^'g^'^) g.xa+ 



dF 



i 9F 
dUiR 



+ 



daiR 



with 



j=0 



(42) 



(43) 



G'"', as above, is the Einstein tensor so, considering the r.h.s. of Eqs.(42,) all the extra 

contributions can be dealt as stress-energy tensor terms. Eqs. (42) arc differential equations 
of order (2 A: H- 4). The stress-energy tensor T^^ contains the kinetic part of the scalar field 
and of the ordinary matter 



(44) 
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The (eventual) contribution of a potential V{(l)) is contained in the definition of F. Prom 

now on, a capital F will indicate a Lagrangian density containing also the contribution of a 
potential V{(j)) while /(</>), f{R), or f{R, Di?) will indicate functions of such fields without 
potential terms. 

By varying the action <S with respect to the scalar field 0, one obtains the Klein-Gordon 
equation 

dF 

.□0=--. (45) 

Dynamics is completely assigned taking into account Eqs.(42), (45) and Bianchi identity 
for Tfj,^. 

Several approaches can be used to treat such equations. For example, by a conformal 
transformation, it is possible to reduce an extended theory of gravity to a (multi) scalar- 
tensor theory of gravity [46,48,52,66]. The theories studied below are specific examples 
coming from the general scheme which we have outlined here. 



C. Scalar— tensor gravity 

The scalar-tensor theories of gravity arc theories, in which both the metric tensor gi^,^ 
and a fundamental scalar field (j) are involved [42]. Their action can be recovered from the 
above (41) with the choice 

F{(j),R) ^ f{(j))R-V{(j)) and e = -1 (46) 

where the coupling f{4>) and the potential V{(p) arc generic functions of the scalar field 
(/). The increasing interest in such theories is mainly due to the success of the inflationary 
paradigm in solving the Standard Cosmological Model shortcomings. Moreover, the presence 
of a scalar field is relevant also in multidimensional gravity, especially in Kaluza-Klein 
theories and in the effective action of String Theory. 

In this framework, the strength of gravity, given by the local value of the gravitational 
coupling, is different from place to place and time to time. For instance, the Brans-Dicke 
theory (the most famous scalar-tensor theory [3]), includes the hypothesis suggested by 
Dirac of the variation of the gravitational coupling with time [43] . 

All these scalar-tensor theories do not satisfy the SEP as the variation of the gravitational 
constant Geff - which is, in general, different from Gn, the standard Newton gravitational 
constant - implies that local gravitational physics depends on the scalar field. Theories 
which have such a feature are called non-minimally coupled theories. 

Some Lagrangians corresponding to particular scalar-tensor theories of gravity are: 

• the effective string-dilaton Lagrangian [40] 

£ = V^e-^(i? + 5''>;^<^;.-A) (47) 

which we shall discuss in details below; 

1 6'^ 

• the Callan-Coleman-Jackiw Lagrangian [67] in which /(0) ^ 2 ~ 12 ' 

• the conformally coupled theory [68] 

C = V^(-^R+^ <?'^>;m'/';. - (48) 
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where V{(p) is generic and f{<p) = — —. 

These Lagrangians arc arc all particular cases of (46). 

In these theories, the gravitational coupling is determined by the form of /(</>). We can 
have two physically interesting situations which could be tested by experiments: 

1) Geff{4>{t))t^oo — * Gn, the standard Newton gravitational constant is asymptotically 
recovered and GR restored; 

2) the possibility that gravitational coupling is not asymptotically constant, i.e. Geff is 
always varying with the epoch and Geff /Geff\now 7^ 0. 

The variability of G can be actually tested by three classes of experiments: 

• the information coming from the Sun dynamics. 

• The observations based on LLR combined with the information obtained studying the 
Sun. As we saw above, LLR consists in measuring the round-trip travel time and thus 
the distance between a transmitter and a reflector located on the Moon; the change 
of round-trip time contains information on the Earth-Moon system. The round-trip 
travel time has been investigated for many years; combining these data with those 
coming from the Sun evolution and the Earth-Mars radar ranging, the more update 
estimation of G/G ranges from 0.4 x 10"-^^ to 10~^^ per year [69,70]. 

• The binary pulsar systems. In order to extract data from his type of system, it has 

been necessary to extend the post-Newtonian approximation, which can be applied 
only in the presence of a weakly gravitationally interacting n-body system, to strong 
gravitationally interacting systems. The estimation of G/G is 2 x 10~^^ per year 
[70,71]. 

Moreover, gravitational Icnsing observations are going to become a further class of tests for 
the variability of Gn [15] 

There is another fundamental issue, coming from the basis of metric theories, to take 
into account scalar-tensor gravity. According to the Mach principle, the local inertial frame 
is determined by the average of the motions of distant astronomical objects [41]. Trying 
to incorporate this principle into metric theories of gravity, Brans and Dicke [3] have con- 
structed a theory of gravity where the total matter distribution influences the strength (and 
the coupling) of gravitational field in each point of spacetime. Therefore, a scalar field <j) 
is introduced as a new gravitational variable together with the standard metric tensor g^^. 
The Brans-Dicke theory is the first and simplest example of a scalar-tensor theory of gravity 
described by the Lagrangian density 

C = ^ (^<PR-uj^^+Cm^ (49) 

where a; is a dimensionless parameter and Cm, as usual, is the Lagrangian density of matter 
including all the non-gravitational fields. 

Let us now perform a conformal transformation g^i, gf^i, = Xg^n, on the Lagrangian 
(49). The field A is proportional to the scalar field (j), i.e. A = Go0 where Go an arbitrary 
constant. The Lagrangian (49) is transformed as 

C = ^ (R + GoCm + GoCx) (50) 

where 
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= - (ot) 

and Cm is the conformally transformed Lagrangian density of matter. Einstein equations 
can be written as 

R^v—^dtivR = Gof^v (52) 

where the stress-energy tensor is the sum of two contributions 

V = T^ir^+A^.(A^-"), (53) 

so that all the contributions due to the scalar field A (or (p) result as a further term in the 
stress-energy tensor. 

This approach constitutes an alternative to GR and several agreements with experimental 

results demonstrated that GR is not the unique possible theory of gravity. 

Moreover, the transformed form of scalar-tensor theory of gravity has the advantage to be 
more manageable, being similar to Einstein standard description. Nevertheless, this form 
presents some drawbacks, e.g. the rest mass is not constant in the conformal transformed 
frame, so the photons trajectories are still geodesies but the massive particles equations 
of motion are modified. This new approach to gravitation has increased the interest in 
scalar-tensor theories of gravity and, especially, in those in which the gravitational constant 
is variable. The Newtonian limit of induced-gravity theories, where the scalar field is non- 
minimally coupled to the geometry, strictly depends on the parameters of coupling and 
self- interaction potential as it has to be by a straightforward PPN parametrization [8] . 

Considering the choice (46), the most general action of a theory of gravity where a scalar- 
field is non-minimally coupled to the geometry is, in four dimensions [73], 



-J 



(54) 



It is possible to obtain the field equations by varying the action (54) with respect to the 
metric tensor field Qixv Therefore, Einstein equations are 

The stress-energy tensor is defined as 

= Tlt^+Tl^^ = -2mTlfJf^ (56) 

i.e. T^v includes the scalar-field and the standard matter terms; in Eq. (55) and Eq. (56), 

TjxJ^'' is the effective stress- energy tensor which contains terms of non-minimal coupling, 
kinetic terms, the potential of scalar field and the usual stress-energy tensor of matter, 
T^™^ , calculated varying with respect to the metric. It is 

+ /(<^);^.-47rGTW} , (57) 

G is a dimensional, strictly positive, constant. Terms coming from the coupling f{4>), 
which were outside T^^^, in Eq. (42), have been assembled. Here Q, defined in Eq. (43), is 



22 



given hy Q = f{^). The standard Newton gravitational constant is replaced by the effective 
coupling 



which is, in general, different from Gat; f{4>) is the generic function describing the scalar field 
coupling. In units with c = 1, Einstein gravity is obtained when the scalar field coupling /((/)) 
is a constant which value is — ^ and G reduces then to the Newton gravitational constant 
Gn [42]. 

The corresponding Klein-Gordon equation is 

D<p-Rf'{<p) + V'{(t>) = 0. (59) 

The derivation of such an equation from the contracted Bianchi identity for r^,y is dis- 
cussed in [42,72,73]. As a general feature, the models described by Eq. (54) can be singularity 
free [74] ; then, there are no restrictions on the interval of time on which the scale factor a{t) 
and the scalar field (f>{t) can be defined. 

Field equations derived from the action (54) can be recast in a Brans-Dicke equivalent 
form by choosing 



^ = f{ct>) u; (<p) = W{^)=V ■ (60) 



The action (54) now reads 



S = j d'^xy/^ 



(61) 



which is nothing else but the extension of the original Brans-Dicke proposal (where w is a 
constant), with ui = uj{ip) plus a potential term W{(p). The sign of u; is chosen according 
to those of coupling and potential. By varying the new action with respect to g^v and the 
new scalar field tp, one obtains again the field equations. The effective stress-energy tensor 
in Eq. (57), in the same units as before, becomes 

Tieff) = _i^rW-^f(^ (P 0,] + 

+ Qi^u^f)- j^9^^i^W{ip) (62) 

and the Klein-Gordon equation (59) becomes 

2 oj {<p) a<p - g"^v-M - g'^'^^-M - + {<fi) = . (63) 

If dip 

This last equation is usually rewritten eliminating the scalar curvature term R with the 
help of Eq. (62), so that one obtains 



(64) 



The minus sign in the denominator comes from the sign convention chosen in the action 
(61). The linearized equations can be derived from the action (54) or, equivalently, from 
(61) as we shall see in the next section. 
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D. Higher-order gravity 



Higher-order theories can be reduced to minimally coupled scalar-tensor ones, and vice- 
versa, by a conformal transformation [45] . The simplest case to be considered is 



C 



-9f{R) 



(65) 



where g is the determinant of metric g^i, and R is the curvature scalar. This Lagrangian is 
recovered from the extended action (41) with the choice F = f{R) and s = 0, that is 



S = Jd^x^ [f{R)+C,n] . 



(66) 



This action depends on a generic function of the Ricci scalar so the Hilbert-Einstein GR 

action is just the linear case sec Eq. (5). Once more, geometric imits 87rG = c= 7i = l 
are used. From the Lagrangian (65) one gets the fourth-order equations 



(67) 



which are fourth-order equations, because of the term /'{Ry^". The prime indicates the 
derivative with respect to R- the standard Einstein equations are immediately recovered if 
f{R) = R. 

The introduction of a new set of variables 



P = f'{R) = f'{9t^i^,dagt^u,dapg^,) 



and 
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P9tJ.i' 



(68) 



is interpreted as the Einstein frame {p, g^v variables), where the word "frame" means 
"change in variables" and the word "Einstein" means that using the transformation 
g (p, g) the non-standard equations (67) are transformed into something closer to the 
standard form, also called the Einstein form. Indeed, the {g^v,p) Einstein equations, with- 
out considering matter terms, are 



G - 1 



lg^u9"^P,aP,0 + 1 gAfiR) - Rf{R)) 



(69) 



Defining 



Inp 



(70) 



one gets 



Gfiv — 



- 2 9ixv(p,a<l 



+ 9i.^V{(l}) 



(71) 



where V(^) is defined as 



V{ct>) 



1 [/(i?)-i?/'(i?)] 
2 



f'\R) 



(72) 



and where the function R = R{p{^)) can be defined thanks to the condition f'{R) ^ 0. 
Eq. (71) corresponds to the Lagrangian 



jC = 



R+7ir''<p,^r -v{ci^) 



(73) 
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which describes the gravitational interaction with a self-gravitating, "minimally" coupled, 
massive scalar field. 



Eq. (67) can also be written in the above Einstein form G^i, = T^^ by defining 

V = . (74) 

In the framework of higher-order-scalar-tensor theories, another choice is extensively 
studied [45,54] 

F = F{R, (j)) any e and = (75) 

with the action 



-9 



F{R,<^)--g^^''ct>,^4>,, . (76) 



The Einstein equations are 

1 

Q V '2 



G^. = ^ (t"" + \{F-gR) + [g.^, - g^^ag] | (77) 



where Q and T^iy are 

OF £ 1 

5 = — and V = 2 ^^■'t^'^-^- - 2 '^^"'^•"l ■ 

The (eventual) contribution of a potential V{<j)) is contained in the definition of F. By 
varying S with respect to the scalar field one obtains the Klein-Gordon equation of the 

form (45). 

Finally, a pure higher than fourth-order gravity theory is recovered with the choice 

F = F{R,nR) e = and = (79) 

which is, in general, an eighth-order theory. If F depends only linearly on OR, one has a 
sixth-order theory. In this model, the action (41) becomes 

d'^x^/^F{R,UR) . (80) 

The Einstein field equations are now 
1 

(81) 

where 

^ dF OF , ^ 

Another generalization of the Hilbert-Einstein action is obtained by including scalars of 
more than second order in derivatives of the metric [75]. The action for these theories has 
the form 



G^'•' = — {g^"'[F - gR] + 2g'i^" - 2g>"' ug - gi^"" [T-^R^ + TUR] + [r^i^R" + T'^Ri^]} 



S = j dTx [f{<i))R + ai iJ' + a2 iJ^'i?^. + as g"" V^,RV.R + 



(83) 



where the a coefficients arc coupling constants and the final suspension points represent 
every other scalar invariants which one can assemble from the curvature tensor, its con- 
tractions and its derivatives. This kind of actions are particularly interesting in one- loop 
quantum gravity since well represent self-gravity interactions and scalar field interactions 
[7]. 
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IV. NEWTONIAN LIMIT IN EXTENDED THEORIES OF GRAVITY 



As discussed in Sec.II, several relativistic theories do not match perfectly GR theoretical 
results but generalize them by introducing corrections to the Newtonian potential which 
could have interesting physical consequences. For example, some theories give rise to terms 
capable of explaining the flat rotation curve of galaxies without using dark matter - for 
instance the fourth-order conformal theory proposed by Mannheim et al. [11]. Others use 
Yukawa corrections to the Newton potential for the same purpose [10]. Besides, indica- 
tions of an apparent, anomalous, long-range acceleration revealed from the data analysis 
of Pioneer 10/11, Galileo, and Ulysses spacecraft trajectories could be framed in a general 
theoretical scheme by taking into account Yukawa-like or higher-order corrections to the 
Newtonian potential [12]. 

In general, any relativistic theory of gravitation can yield corrections to the Newton po- 
tential which, in the PPN formalism, could furnish tests for this theory [8]. In this section, 
the Newtonian limit is investigated for fourth-order and scalar-tensor theories showing that 
several new features come out with respect to the corresponding GR Newtonian limit. 

A. Recovering the Newtonian limit in fourth-order theories of gravity 

Particularly interesting, among the fourth-order theories of gravity, are the so called 
i?^-theories [49] , where the inclusion of terms proportional to iJ^^i?'*'' and in the gravi- 
tational action gives rise to a class of effectively multi-mass models of gravity [9] ; such terms 
also produce a stabilization of the divergence structure of gravity, allowing it to be renor- 
malized through its matter couplings. The empty space solutions of Einstein equations, 
J^^ll' _ .^-^^ g^jgQ solutions of the field equations derived from actions like / \/— g R^i.R^'^ 
and / ^/^B?, but the classical tests of GR are not automatically satisfied. In fact, al- 
though e.g. the Schwarzschild solution satisfies the empty space equations, it is not the one 
that couples to a positive definitive matter distribution; those solutions of purely fourth- 
derivative models which do couple to a positive matter source are not asymptotically flat at 
infinity. 

Let us consider only models derived from actions including both the Hilbert action 



/ ^/—g R and the four derivative terms. There are only two independent possible additions, 

because of the Gauss-Bonnet relation in four-dimensions, so there is only a two-parameter 
family of field equations without cosmological term (see also [38]). 
The higher-order action is, in this case, 



where fc^ = 327rGAr, and a, (5 and 7 are dimensionless numbers. It turns out that the 
correct physical value for 7 is 2, as in the Einstein theory, restoring the standard coupling 



As we shall see, the static spherically symmetric solutions to the field equations derived 
from Eq. (84) either reduce asymptotically to the sum of a Newtonian and two Yukawa po- 
tentials, or they are unbounded at infinity and must be eliminated by some proper boundary 
conditions. The masses, in these Yukawa potentials, are only weakly constrained by the ob- 
servational evidence to be > 10^^ cwr^ (in geometrical units). Although the magnitude of 
these effects is negligible at laboratory distances, there are some interesting features which 
could result interesting at astronomical distances beyond the Solar System. Coupling the 
linearized theory to a pressurized fluid distribution shows that the coefficients of the Yukawa 
potentials depend on the pressure and the size of the distribution. This shows that Birkhoff 's 





(84) 



(167rGjv) 



-1 
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theorem is not valid in these models. 

The Yukawa potentials and the failure of BirkhofF theorem in the static case arc comple- 
mented by the existence of " massive" radiation terms in the dynamical field. The analysis 
of the field dynamics, both in the action and in the interaction between sources, shows 
that the gravitational field has eight degrees of freedom, two corresponding to the mass- 
less spin-two graviton, five to a massive spin-two excitation and the last one to a massive 
scalar. This decomposition is carried out into various individual spin components using the 
TT decomposition of a symmetric tensor. The massive excitations are separated from the 
massless ones by introducing auxiliary oscillator variables, in terms of which the Lagrangian 
may be separated into a sum of Lagrangians for massive scalar-fields. The massive spin-two 
field comes with a minus sign relative to the other fields, both in the oscillator variable La- 
grangian and in the radiation. Classically, the corresponding excitation has negative energy; 
therefore, there is a breakdown of causality which can cause superpositions of waves with 
group velocities greater than the speed of light. 

The field equations obtained by the action (84), including matter, are 

H^, = (a - 2/3) R,^,, - a - Q a - 2 iJj^,.^ + 2 a RP^R^^^x (85) 

-1iiRR^,,-\^ g^, {a R'^^R^x - /3 i?') + 7 R^^-^l R ^ ^ (86) 

Eq. (85) and Eq. (86) are linked by generalized Bianchi identities Hf^'^ = 0. In order to 
extract the physical content of Eq. (85) for the static spherically symmetric gravitational 
field, one can work in Schwarzschild coordinates, see Eq. (1). The above equations can be 
linearized assuming A{r) = 1 + V{r) and B{r) = 1 + W{r) and keeping only terms linear in 
V and W. 

The general solution to the linearized field equations (85), i.e. Hj^^, = 0, is 

V = C+ — + C^+- + C^-- + C0+^ + (87) 

y y y 

rilfi „m2 r c—m-2 r „mo r o—ma r 

W = -— -C'+- C^- + C0+ + + 

y rp rp rp 

1 1 p-mar 

+ -C^+TOa e'"^ - -C^-mz e"™^ - C°+too e™" + C^-mo (88) 

2 2 r 

where V is the gravitational (Newton-like) potential and m2 and mo are related to the 
coupling constant a, /3, 7 and k as 



^2 = a/-^2, and mo = . L,^^' . , o , (89) 



„fc2' — ""u Y2(3/3-a)fc2 

C, C^'°, C^+, C^~, C°+, C°~ are arbitrary constants. The rising exponentials r~^e"^'^ occur 
in Eq. (87) and Eq. (88) because they are solutions of the static Klein-Gordon equation (V^ — 
m^) (f) = 0. As in the Klein-Gordon case, they must be eliminated by imposing boundary 
conditions at infinity. Counting the number of parameters that must be determined by 
source coupling and boundary conditions at infinity, the parameter C should be discarded, 
since it may be absorbed into a trivial rescaling of the time coordinate t. Consequently, the 
general solution to the linearized static spherically symmetric equations is a five-parameter 
family, before the imposition of boundary conditions. 

In the case of a massive point particle with T^^ = 6° S'lM S^{x), the gravitational potential 
takes the form 

V = -- \-- — . (90 

87r7r 6777 r 24777 
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The comparison at infinity with the Newtonian result V = —2 G M shows that the 

correct physical value of 7 is 2. The higher-derivative terms can produce an appreciable 
effect also at small distances from the source of the field, at a scale set by 1712 and toq. 
At the origin, the Newtonian term 1/r is cancelled and Eq. (90) tends to the finite value 
fc2 M (24 TT 7)"^(mo - 4 777,2). 

If, instead of a point particle, one couples to an extended source with some structure, the 
peculiarities of the static field coupling in these models will be revealed. One can consider 
now a spherical distribution of a perfect fluid with an internal pressure that is maintained 
by an elastic membrane. The stress tensor for this source is 



T = 



(P 

P[l 







\l5{r 





0]^ 






P[l-\l5{r-l)]r'^ sin^ 9 

3M(47rZ3)-iy 



(91) 



where I is the radius of the source, M is its total energy and P is the internal pressure. The 
delta functions in T00 and T^^ correspond to the membrane tension and are necessary to 
ensure equilibrium of the source, given by the Bianchi identities Tl^" = 0, which read, in 
spherical coordinates, as 



{TrrY + 2r-'^Trr-rT0e-r sin'^eT^^, = 



(92) 



The exterior field, for the source Eq. (91), is 



Trext 
^ fluid 



H Si S2 



7r 



27rr 



(93) 



where si and S2 are defined as 



si = 



S2 



M 



I cosh(7772 sinh(7n2 



- P 



sinh(7n2 



ml 



+ 



cosh(r772 I) P sinh(7n,2 /) 



3m2 



- P 



M 
47r/3 

sinh(mo I) 



I cosh(777o /) sinh(7r7o /) 



ml 



ml 



I cosh(mo I) P sinh(mo I) 



mt 



mf, 



+ 



3 mo 



(94) 



In Eq. (93), the existence of the quantities 77t2 and mo with the dimensions of mass has 
enabled the gravitational coupling to involve the source pressure and radius, as well as the 
total energy M. In the limit I 0, the pressure dependence drops out and Eq. (93) tends 
to the potential (90) for a point particle. The three-parameter family of static solutions - 
Eq. (93) - shows that the Birkhoff theorem does not hold for this models; this is the case for 
the Yukawa potentials in Eq. (90). Eq. (93) corresponds to the virtual exchange of massive 
particles and has a real counterpart in massive spin-two and spin-zero radiation. Eq. (90) 
and Eq. (93) give an acceptable Newtonian limit only for real m2 and mo, so that there are 
only falling exponentials and not oscillating 1/r terms at infinity. This corresponds to the 
absence of tachyons (both positive and negative energy) in the dynamical field. It is also 
possible to remove the two massive fields by picking combinations of a and P that make m2 
or mo infinite, choosing, respectively, a = or a = 3/? in Eq. (85). 

The static field (90) may be considered as the gravitational potential of a star (or a generic 
spherical body) and observations could be used to compute bounds on the values of m2 and 
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mo- However, astronomical tests in Solar System are useless here. In fact, by considering 

the motion of Mercury, the orbital precession is known, up to now, to about one part in 10^. 
The corrections due to the higher-derivative terms are of order e~""^, where, for r, one can 
use the radius of Mercury orbit 5 x 10^ km. This gives a lower bound on the masses of 
the order ^ 5 x 10^"'^^ cm^^. Since there is no discontinuity in the coupling to light in the 
limit m2, mo oo, measurements of light bending by the Sun do no better. Laboratory 
experiments on the validity of the Newtonian force seem much more promising, as they 
can seek agreement on the meter scale. Yet, high-precision Cavendish-type experiments 
have to be improved. On the other hand, such results could provide interesting clues at 
extragalactic scales to explain flat-rotation curves of galaxies [11] but the debate is open 
and final tests should be achieved. 

B. Newtonian limit in scalar-tensor gravity 

Also in the case of scalar-tensor gravity, the Newtonian limit gives non-standard gravi- 
tational potentials with exponential contributions directly depending on the parameters of 
scalar field coupling and potential. Before starting with our analysis, we need a choice for 
the up to now arbitrary functions /{(p) and V{(f)) used in action (54). A rather general 
choice is given by 

m = i<r , (95) 

y(</.) = , (96) 

where ^ is a coupling c;onstant, A gives the self interaction potential strength, m and n are 
arbitrary parameters. This a physically motivated choice, in agreement with the existence 
of a Noether symmetry in the action (54) as discussed in [42,73]. Furthermore, several 
scalar -tensor physical theories (e.g. induced gravity) admit such a form for f{(p) and V{(p). 

It is instructive to develop in details the calculations in order to see how extra contri- 
butions come out. In order to recover the Newtonian limit, we write, as above, the metric 
tensor as 

fl/iii. = Vni^ + , (97) 

where ry^jy is the Minkoskwi metric and h^i, is a small correction to it. In the same way, we 
define the scalar field ip as a perturbation, of the same order of the components of /i^y, of 
the original field (j), that is 

<P = ipo + i', (98) 

where (po is a constant of order unit. It is clear that for ipQ = 1 and tjj = Einstein GR is 
fully recovered. 

To write in an appropriate form the Einstein tensor G/^u, we define the auxiliary fields 

hj^i/ = hjj,i/ — T/i^i/h , (99) 

and 

(T„ = hc.f^,jri^'' . (100) 
Given these definitions, to the first order in hf^v, we obtain 
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G;^!/ = {nr/^pi/ + riniy'^a,0V"'^ - - <^iy,ii] , (101) 

where = rj^'^d^d^. We have not fixed the gauge yet. 

We pass now to the right hand side of field equations (55), namely to the effective stress 
energy tensor. Up to the second order in ip, the coupHng function f{^) and the potential 
V{(j)), by using Eqs.(95) and (96), become 

f{^) + m^^-i V + ^K r^ , (102) 



n 



(n-1) 



y(<^) :^X(^^+n ^^ V + ^ ^ ' ^oV ] ■ (103) 



To the first order, the effective stress-energy tensor becomes 

m+n 

io _ 

2^ ^^'^ ^ 



T^, = -mvp^—i rj^^D^^ + mipl"'-'^p,^, - t?;.. - (47rG)^ T^, , (104) 



and then the field equations are 



X m+n m 

^r,,. + (4.G)^ 



We can eliminate the term proportional to by choosing an appropriate gauge. In fact, 
by writing the auxiliary field cfa, given by Eq.(lOO), as 

C^a = m^^"-V,a, (106) 

field equations (105) read 

□^V-'^Vo""V°r,V'- ^^^??M. + (8 7rG)^ V (107) 

By defining the auxiliary field with components /i^j/ as 

h^,u = hfiu - m(fl"'''^T]ij,utjj , (108) 
the field equations take the simpler form 

□r,V = 'y^*- + (SttG)^ T^. (109) 

The original perturbation field h^v can be written in terms of the new field as (with h = 

h,,u = h^^-^r]^Ji-rnipl"'~'^r}^^'il) . (110) 

We turn now to the Klein-Gordon Eq.(59). can be written, from the linearized Klein- 
Gordon equation, in terms of the matter stress energy tensor and of the potential term. If 
we calculate the scalar invariant of curvature R = g^^^R^v from Eq.(57), we find 
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— '^^ 9"^<l>;o.<t>;0 - 2 V{cf) - 3 □/(</.) - 47rGT ) + V'{<t>) = , 



/(<A) 

and, to the first order, it reads 

^ A(n-2m)(n-l)yg-^ ^ A(2m - n)./pg-^ 



A-kGw? 



1 - 3emVo 



1 - 3 ^mV""' (1 - 3 emV™"')<Po 



(Ill) 



(112) 



We work in the weak- field and slow motion limits, namely we assume that the matter stress- 
energy tensor T^^, is dominated by the mass density term and we neglect time derivatives 
with respect to the space derivatives, so that — > —A, where A is the ordinary Laplacian 
operator in flat spacetime. The linearized field equations (109) and (112) have, for point-like 
distribution of matter*, which is p{r) = MS{r), the following solutions: 
for n 7^ 2m, n ^ 1, we get 
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cosh(pr) , 

.—pr 



1 — 3 ^m'^Lpf) 



47r(/?o 
n-1 



cosh(pr) , 



where the parameter p is given by 



2 _ A(n-2m)(n- 1)<^^ 



l-3^m2(^S' 



For n = 2m, we obtain 
hoQ ~ 



(4^G)^fr(l-4CmVr')' 
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(113) 

(114) 
(115) 

(116) 

(117) 
(118) 
(119) 



2m- 1 

■f 

values m = 1, n 
potentials, i.e. / ' 



where A = mLpfj^^^^ipQ and -00 are arbitrary integration constants. Let us note that the 
2 and m = 2, n = 4 correspond to the well known couplings and 
,V (fp' and f ^ (jy^ , V ^ (f)'^ , respectively. 



*To be precise, we can define a Schwarzschild mass of the form 

M = j{2TS-Tp^d'x. 
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Finally for n = 1, we obtain 

"(47rG)^^(l-4CmVo""') 



hii ^ Sii 



e(i-3^mvr') 
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(121) 



(122) 



If we demand the (0, 0)-component of the field Eq.(109), when A = 0, to read as the usual 
Poisson equation (that is nothing else but a definition of the mass) 



A$ = AttGnP, 



(123) 



where linked with the metric tensor by the relation /iqo = 2 is the Newtonian potential, 
we have to put 

4^mVjr"^^ 



r ^0 
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G. 



(124) 



2^ Vl-3em2<^S^-^ 

We may now rewrite the nonzero components of hf^„ and the scalar perturbed field. Let us 
take into account, for example, Eqs. (117)-(119). We get 



hoo — — 



2GnM 47rAv5o 
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(125) 



(126) 



(127) 



where the Newton constant explicitly appears. Similar considerations hold in the other 
cases. 

What we have obtained arc solutions of the linearized field equations, starting from the 
action of a scalar field nonminimally coupled to the geometry, and minimally coupled to the 
ordinary matter. Such solutions depend on the parameters which characterize the theory: 
^,771, n, A. The results of Einstein GR are obtained for / = /o, with /o negatively-defined 
due to the sign choice in the action (54). As we can easily see from above, in particular 
from Eqs.(125)-(127), we have the usual Newtonian potential and a sort of cosmological 
term ruled by A which, from the Poisson equation, gives a quadratic contribution. 

We consider now the Brans-Dicke-like action (61) where w = w (t/j). It is actually simple 
to see, from the Eqs. (62) and (63), that, if we want to limit ourselves to the linear approxi- 
mation, we may consider as well lo =constant. The link with the results that follows from 
the action (61) are given by the transformation laws (60), that is 

^(<pW) = ^ (128) 

The potential term W {(f) in the linear approximation, behaves as V {(/)). Results in the 
approaches given by the actions (54) and (61) are completely equivalent. 
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C. An application: Newtonian limit of String-Dilaton Gravity 

The above approach can be applied to the case of string-dilaton gravity which is a remark- 
able example of scalar-tensor theory which have to be investigated also in the Newtonian 
limit. 

Let us shortly outline why such a theory is so interesting from cosmological and astro- 
physical viewpoints. 

String-dilaton gravity seems to yield one of the most promising scenarios in order to solve 
several shortcomings of standard and inflationary cosmology [76]. First of all, it addresses 
the problem of the initial singularity, which is elegantly solved by invoking a maximal space- 
time curvature directly related to the string size [77]. Besides, it introduces a wide family 
of cosmological solutions which comes out thanks to the existence of a peculiar symmetry 
called duality, which holds at string fundamental scales, as well as at cosmological scales 
[40]. In practice, if a{t) is a cosmological solution of a string-dilaton model, also a~^{t) 
has to be one by a time reversal t — » —t. In this case, one can study the evolution of the 
universe toward t +oo as well as toward t —oo. The junction of these two classes of 
solutions, at some maximal value of curvature, considered as branches of a (up to now not 
found) general solutions, has to be in agreement with inflationary paradigm and solve the 
initial values and singularity problems of standard cosmological model [77] . 
String-dilaton cosmological models come out from the low energy limit of (supcr)string the- 
ory which can be considered as one of the most serious attempt to get the great unification of 
all forces in Nature in the last thirty years. This theory avoids the shortcomings of quantum 
field theories due, essentially, to the point like nature of particles the renormalization - 
and includes gravity in the same conceptual scheme as the other fundamental interactions 
- the graviton is just a string mode like the other gauge bosons [78] . 

However, despite the quality and the variety of theoretical results, the possibility to test 
experimentally the full predictions of the theory is not yet available. The main reason for 
this failure is that the Planck scale - ~ 10^^ GeV -, where the string effects become relevant, 
is unreachable in today's high-energy physics experiments. Therefore, cosmology and astro- 
physics remain the only open ways for observational investigations and detecting remnants 
of primordial processes could lead to test this theory. Furthermore, a lot of open questions 
of astrophysics, like dark matter, relic gravitational-wave background, large-scale structure, 
primordial magnetic fields and so on could be explained by strings and their dynamics. 
The key element of string-dilaton gravity [76,79], in low-energy limit, is the fact that a 
dynamics, consistent with duality, can be implemented only by taking into account massless 
modes (zero modes) where the scalar mode (the dilaton) is non-minimally coupled to the 
other fields. The tree-level effective action, in general, contains a second-rank symmetric 
tensor field (the metric), a scalar field (the dilaton) and a second-rank antisymmetric ten- 
sor field (the so-called Kalb-Ramond universal axion). Such an action can be recast as a 
scalar-tensor theory, e.g. induced gravity, where the gravitational coupling is a function 
of the dilaton field [72,73,76]. Then, it is legitimate to study the Newtonian limit of the 
string-dilaton gravity to see what is its behavior in the weak-field and slow-motion approx- 
imations. This approach is useful to investigate how string-dilaton dynamics could affect 
scales shorter than cosmological ones. The issue is to search for effects of the coupling and 
the self-interaction potential of the dilaton at scales of the order of Solar System or Galactic 
size. 

Calculations can be developed in the string frame since one wants to understand the role 
of a dilaton-non- minimal coupling in the recovering of the Newtonian limit. In the string- 
dilaton theory, it is necessary to introduce a scalar potential with a power law expression of 
the form V{(j)) = X (sec [80] for details). 

The tree-level string-dilaton effective action, the lowest order in loop expansion, contain- 
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ing all the massless modes, without higher-order curvature corrections {i.e. without the 
Gauss-Bonnet invariant) is 



S = - 




(129) 



where R is the Ricci scalar, </> the dilaton field, V{(j)) the dilaton self-interaction potential, 
-ff/ji/a = c^[/i-Bi/a] is the full antisymmetric derivative of the Kalb-Ramond axion tensor and 
Cm is the Lagrangian density of other generic matter sources. The theory is formulated in 
(d+l)-dimensions, d being the number of spatial dimensions, and Ag the string fundamental 
minimal length parameter. To the lowest order, the effective gravitational coupling is given 
by 



As e 



0/2 



(130) 



where is the Newton constant and Ap is the Planck length. Units have been chosen such 
that 2Ag~^ = 1 so that e'^ is the {d + l)-dimensional gravitational coupling. 

Clearly, the string- dilaton effective action (129), apart the Kalb-Ramond term, is a par- 
ticular case of the general scalar-tensor action (54) where the transformations 



<p — > exp[-(/)] , /(0) — > i exp[- 



y(^) _^v^(0)exp[ 



(131) 



are performed. 

The field equations are derived by varying the action (129) with respect to g^^, 4>, and 
B^i,. One gets, respectively. 



1 

2 SiJ.i' 



R + 2V^4>- {V4>f +V-V' -— H^,„ 



(132) 

(133) 
(134) 



where, as above, G^v is the Einstein tensor, T^i, the stress-energy tensor of matter sources 
and V = dV/d(j). 

Let us assume that standard matter is a perfect fliiici minimally coupled to the dilaton. 
Otherwise, in the non-minimally coupled case, a further term should be taken into account 
in Eq. (133). The above equations form a system of tensor equations in {d + 1) dimensions 
which assigns the dynamics of d, and B^i,. 

Now let us consider the weak field approximation in order to derive the PPN limit of 
the theory. Obviously, all the invariants of the full theory are not preserved if equations 
are linearized. For example, duality is lost in the linearized solutions. However, this is not 
a problem in the present context since there is the assumption of a regime well-separated 
from early singularity where duality is adopted to solve cosmological shortcomings. Here, 
the aim is to check if rem,nants of primordial string-dilaton dynamics are detectable at our 
nearest scales (Solar System or Galaxy). 

Moreover, the weak field approximation for the axion gives only second order terms with 
respect to ft-^j/ and ip in the field equations, so that its contribution in the following consider- 
ations can be discarded. A physical interpretation of this fact could be that the production 
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of primordial magnetic fields, considered as "seeds" for the today's observed large magnetic 
fields of galaxies [81] is a second order effect if it is due to Hf^ua- 

Using the approximated expressions (97) and (98) for the fields and the natural assump- 
tion V{(j)) = X (j/^ for the scalar field potential, the field equations (132) and (133) become 
respectively 

^■qhni, + T]^^ <T„ " - (<7^,^ + (T^,^) - 2 V, ^l/ + 2 T]^^ D^V + A {V^u + h^,,) (^q + 

+Anv'^^?7^,V - -e'^°T^. (135) 
2U^xl,+ ]^U^h + aJ^ + X^p^[nipQ^ -n{n-l) ^q^] + A [1 - n ^] = • (136) 

The field equation (134) has been discarded since it gives only terms higher than the 
linear order. The term proportional to V.jul/ can be eliminated by choosing an appropriate 
gauge. In fact, by writing the auxiliary field <Ja as 

= (137) 

Eq. (135) reads 

~ -e^° T^, . (138) 

In this approximation, the terms in h^v and in the brackets can be neglected as /i^ji^ <C 
'qnv and ij) ^1. Eq. (138) then becomes 

^rthf.u+Vt.uU^i^ + XLp'^ri^^ ~ -e^oT^^. (139) 

By defining the further auxiliary field 

h^^ = hi^^+rji^^ip (140) 

one gets the final form 

□r,V + ^^^0 -e"'^°V • (141) 

The original perturbation field h^i, can be written in terms of the new field, i.e. 

hij,u = h^u-^ Vfj-i' h + r]ij,^tp (142) 

with h = 'q^'^h^y. 

Let us turn now to the Klein-Gordon equation (136). With some algebra, it can be recast 
in the form 

{u^ + c\) V ^ ^T-$o (143) 

where T is the trace of the stress-energy tensor of standard matter and the constants are 

cl = X^^ [n{n-l)^f -n^^^] , *o = {^-n^^^) (144) 

The linearized field equations (141) and (143) have, for point-like distribution of matter, 
the following solutions 
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hoo ~ ^^^^{1 ~ e-^i n + C2 r2 + C3 cosh(cir) (145) 

r 

hii ~ -i^^(l + e-^i '■) - C2 - C3 cosh(cir) (146) 
r 

^ — — + C3 cosh(ci r) (147) 

where the coefficients C2 and C3 are 

C2 = 27rA(/p^ C3 = . _i (148) 

Here the functions (145) and (146) are the non-zero components of h^v while Eq. (147) 

is the perturbation of the dilaton. Standard units have been restored. 

When A = and in the slow motion limit, the (0, 0)-component of the field Eq. (141) 
reduces to the Poisson equation (123). 

Solutions of the linearized field equations have been obtained, starting from the action 
of a scalar field non-minimally coupled to the geometry, and minimally coupled to the 
ordinary matter (specifically the string-dilaton gravity). Such solutions explicitly depend 
on the parameters (fQ,n,\ which assign the model in the class of actions (129). At this 
point, string-dilaton gravity is completely equivalent to any scalar-tensor theory of gravity 
and the physical considerations are the same. 



V. DISCUSSION AND CONCLUSIONS 



As a general remark, the Newtonian limit of ETG strictly depends on the parameters 
of the scalar-field couplings, the self-interaction potentials and the higher-order coupling 
terms as it has to be by a straightforward PPN parametrization [8] . 

As we have seen, the self interactions (scalar or gravitational) give rise, in any case, to 
corrections to the Newtonian potential. Such corrections are, in general, constant, quadratic 
or Yukawa like terms. Due to these facts, ETG could be the potential natural candidates 
to solve several astrophysical problems as the flat rotation curves of spiral galaxies [9,11,82] 
(and then the up to now not revealed presence of consistent bulks of dark matter) , the ap- 
parent, anomalous long-range acceleration revealed by several spacecrafts [12] . Essentially 
the corrections depends on the strength of the couplings and the "masses" of the scalar fields 
(where, also higher-order curvature terms can be considered under the standad of "scalar 
fields" using conformal transformations). It is worth noting that we get always scale lengths 
where Einstein GR can be, at least approximately, recovered. This fact could account why 
measurements inside the Solar System confirm such theory while outside of it there are 
probable deviations [12]. 

We can make these considerations more precise discussing some of the PPN solutions 
which we found above. For example, let us take into account the Newtonian-like potential, 
derived in PPN approximation, from string-dilaton gravity. 

As shown, the solution given in Eq. (145) can be read as a Newtonian potential with 
exponential and quadratic corrections, i.e. 

^ _^JvM _ ^ ^ ^2 ^ ^ cosh{ci r) . (149) 

r 2 2 

In general, most of ETG have a weak field limit of this form (see also [52,9,82]), that is 

GnM 

— i -I- 



$(r) 



(150) 
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where Gjv is the value of the gravitational constant as measured at infinity and rk is the 
interaction length of the k-th component of the non-Newtonian corrections. The amplitude 
ak of each component is normalized to the standard Newtonian term and the signs of the ak 
coefficients indicate if the corrections are attractive or repulsive [8]. Besides, the variation 
of the gravitational coupling is included. As an example, let us take into account only the 
first term of the series in Eq. (150) which is usually considered as the leading term (this 
choice is not sufiicient if other corrections are needed). One has 

^r) = -^E^i [l + aie-'^/'-i] . (151) 
r L J 

The efi'ect of non-Newtonian term can be parameterized by (ai, ri). For large distances, 
at which r ^ ri, the exponential term vanishes and the gravitational coupling is Gat. If 
r <C ri, the exponential becomes of order unity; by differentiating Eq. (151) and comparing 
with the gravitational force measured in laboratory, one gets 



Glab = G 



N 



1 + ai ( 1 + — ) e 



^ -r/ri 



Gjv(l+ai) (152) 



where Giab = 6.67 x 10~^^ Kg~^m^s~^ is the usual Newton constant measured by the 

most recent Cavendish-likc experiments. Of course. Gat and Giab coincide in the standard 
gravity. It is worthwhile noticing that, asymptotically, the inverse square law holds but 
that the measured coupling constant differs by a factor (1 + ai). In general, any correction 
introduces a characteristic length that acts at a certain scale for the self- gravitating systems. 
The range of of the fcth-component of non-Newtonian force can be identified with the 
mass mfc of a pseudo-particle whose Compton's length is 

rk = ^ . (153) 
rrifec 

This means that, in the weak energy limit, fundamental theories which attempt to unify 
gravity with the other forces introduce, in addition to the massless graviton, particles with 
mass which carry the gravitational force [83] . 

There have been several attempts to constrain rk and ak (and then mk) by experiments 
on scales in the range 1cm < r < 1000 km, using totally different techniques [84-86]. The 
expected masses for particles which should carry the additional gravitational force are in 
the range 10~^^eV < ruk < 10~^eV. The general outcome of these experiments, even 
retaining only the term fc = 1 , is that a "geophysical window" between the laboratory and 
the astronomical scales has to be taken into account. In fact, the range 

[ail ~ 10-2, n ~ lO^^lO^m (154) 

is not excluded at all in this window by current measurements. An interesting suggestion has 
been made by Fujii [87] , which proposed that the exponential deviation from the Newtonian 
standard potential (the "fifth force") could arise from the microscopic interaction which 
couples to nuclear isospin and baryon number. 

The astrophysical counterpart of these non-Newtonian corrections seemed ruled out as 
recent experimental tests of GR "exactly" predict the Newtonian potential in the weak 
energy limit, inside the Solar System. Besides, indications of an anomalous, long-range 
acceleration have been found during the data analysis of Pioneer 10/11, Galileo, and Ulysses 
spacecrafts (which are now almost outside the Solar System) and made these Yukawa-like 
corrections come again into play [12]. Furthermore, Sanders [10] reproduced the flat rotation 
curves of spiral galaxies by using 
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ai = -0.92, ri~40kpc. (155) 

His main hypothesis is that the additional gravitational interaction is carried by an ultra- 
soft boson whose range of mass is rrii ~ W~^^ 10~^^eV. The action of this boson becomes 
efficient at galactic scales without the request of enormous amounts of dark matter to sta- 
bilize the systems. In addition, Eckhardt [86] uses a combination of two exponential terms 
and gives a detailed explanation of the kinematics of galaxies and galax;y clusters, without 
dark matter models, using arguments similar to those of Sanders. 

It is worthwhile noticing that both the spacecrafts and galactic rotation curves indications 
are outside the usual Solar System boundaries used to test GR. However, the above authors 
do not start from any fundamental theory in order to explain the outcome of Yukawa cor- 
rections. In their analysis, Yukawa terms are purely phenomenological. 

Another important remark in this direction deserves the fact that some authors 
[88] interpret the recent data coming from cosmic microwave background experiments, 
BOOMERANG [89] and WMAP [59], in the framework of modified Newtonian dynam- 
ics, without invoking any dark matter model. 

All these facts point toward the line of thinking that "corrections" to the standard gravity 
have to be seriously taken into account. Let us now turn back to the above solutions (145) - 
(147), in particular to the gravitational potential (149). It comes out from the weak field 
limit (PPN approximation) of a string-dilaton effective action (129). The specific model 
is singled out by the number of spatial dimension d and the form of the self-interaction 
potential V{(j)) where the quite general class V{(t)) = Xcp""' has been assumed. 

Without loosing of generality, one can assume <fo = ^- This means that for = 1 the 
standard gravitational coupling is restored in the action (129). However, the condition 
ip must hold. For the choice n = 3 in Eq.(148), one has 

^r) 2. -^IL^{l-e-^n + ^r' (156) 
r 2 

where, beside the standard Newtonian potential, two correction terms are present. Due to 
the definition of the constants ci,2 their strength directly depends on the coupling A of the 
self-interaction potential V{(p). If A > 0, one finds, from Eq. (144) and Eq. (148), that the 
first correction is a repulsive Yukawa-like term with ai = —1 and ri = c^^ = X^^^"^. The 
second correction is given by a kind of positive-defined cosmological constant C2 which acts 
as a repulsive force^ proportional to r. On the other hand, if A < 0, the first correction is 
oscillatory while the second is attractive. 

From the astrophysical point of view, the first situation is more interesting. Assuming that 
the dilaton is an ultra-soft boson which carries the scalar mode of gravitational field, one gets 
that the length scale ^ 10^^-^ 10^^ cm, is obtained if its mass range is m ^ 10~^^-^ 10~^*cV. 
This length scale is necessary to explain the flat rotation curves of spiral galaxy. The second 
correction to the Newtonian potential can contribute to stabilize the system being repulsive 
and acting as a constant density a sort of cosmological constant at galactic scales. In 
general, if ai ~ —1 the flat rotation curves of galaxies can be reconstructed [10]. If the 
mass of the dilaton is in the range 10~^^eV < m < 10~^eV the "geophysical windows" 
could also be of interest for experiments. Finally the mass range m ^ 10~^ 10~"'^'^eV 
could be interesting at Solar System scales (for the allowed mass windows in cosmology see 



^If in the Poisson equation there is a positive constant density, it gives rise to a repulsive quadratic 
potential in r and then to a linear force. A positive constant density can be easily interpreted as a 
cosmological constant. 
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Ref. [90]). Similar analysis can be performed also for other values of n corresponding to 
different models of the class defined in Eq. (129). 

In conclusion, the weak energy limit of string-dilaton gravity has been derived, showing 
that in this case the Newtonian gravitational potential is corrected by exponential and 
quadratic terms. These terms introduce natural length scales which can be connected to 
the mass of the dilaton. If the dilaton is an ultra-soft boson, one can expect observable 
effects at astrophysical scales. If it is more massive, the effects could be interesting at 
geophysical or microscopic scales. 

This is just an example, but a similar analysis can be developed also for other ETG and 
the clues are that wide classes of gravitational phenomena, at microscopic, macroscopic and 
astrophysical scales, could be interpreted by enlarging the scheme of GR. 
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